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We continue the investigations of the 16 O ground state using the coupled-cluster expansion 
[cxp(S)] method with realistic nuclear interaction. In this stage of the project, we take 
into account the three nucleon interaction, and examine in some detail the definition of 
the internal Hamiltonian, thus trying to correct for the center-of-mass motion. We show 
that this may result in a better separation of the internal and center-of-mass degrees of 
freedom in the many-body nuclear wave function. The resulting ground state wave function 
is used to calculate the "theoretical" charge form factor and charge density. Using the 
"theoretical" charge density, we generate the charge form factor in the DWBA picture, 
which is then compared with the available experimental data. The longitudinal response 
function in inclusive electron scattering for 16 is also computed. 

PACS numbers: 21.60.Gx,25.30.-c,21.10.Ft,27.20.+n 
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I. INTRODUCTION 

The coupled-cluster expansion, also called the exp(S) 
method, was developed about 40 years ago by Coester [0. 
It was then applied to finite nuclei quite successfully 
by Kummel and the Bochum group || using a repre- 
sentation of the wave function in coordinate space to- 
gether with common interactions of that time. While 
the method in general is viewed as essentially exact, ap- 
proximations are introduced stemming from truncations 
in the coupled cluster equations, as well as truncations 
in the model space. The results for the binding energy 
were well above the Coester-line and were taken as evi- 
dence for the presence of three-nucleon and higher order 
interactions. 

While few further developments took place since that 
time, the ongoing expansion in computer power as well 
as the availability of more sophisticated nucleon-nucleon 
interactions suggests that one should go back to these 
precise methods. During the past several years we have 
reexamined the coupled cluster approach and applied it 
to the spherical nucleus 16 using the Argonne vl8 po- 
tential |3| . We combine the mean field approach together 
with the coupled cluster expansion in a harmonic oscil- 
lator basis of 50 Tioj (which constitutes our truncated 
Hilbert space) and going up to the level of ApAh clusters. 

Further, we have expanded the formulation to include 
the Urbana IX three-nucleon interaction. The results 
are in reasonable agreement with the experiment and 



make the subject of the present paper. These results 
obtained using the coupled cluster expansion should be 
directly comparable with those obtained using the Vari- 
ational Monte Carlo (VMC) j| or the Green's Function 
Monte Carlo (GFMC) methods, using the same interac- 
tion. Especially the GFMC method with the Argonne 
v 18 potential and the Urbana IX three-nucleon interac- 
tion has proved quite accurate in nuclei with A < 8 [0, 
even though a certain uncertainty in the three-nucleon 
interaction model is responsible for a systematic lack of 
binding in the heavier systems. These difficulties will 
be addressed however in the future, by considering more 
sophisticated three-body interaction models. 

In section [ll] we briefly review the coupled-cluster ap- 
proach as outlined in when the Hamiltonian includes 
only the two-body interaction. We continue this discus- 
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sion in section [II with comments pertaining to the in- 
clusion of the center-of-mass Hamiltonian, which is sup- 
posed to constrain the center-of-mass degrees of motion 
in the ground state wave function. In section ^ we out- 
line the changes necessary in order to include the three- 
nucleon interaction via a density-dependent approxima- 
tion. In section [v| we show various observables calcu- 
lated in the ground state for 16 O. Finally, in section VI 
we present our conclusions and a brief outlook. 



II. COUPLED-CLUSTER METHOD: A REVIEW 

For a spherically symmetric nuclear system, consist- 
ing of both protons and neutrons, the total Hamiltonian 
is given as the sum of a nonrelativistic one-body kinetic 
energy, a two-nucleon potential and a three-nucleon po- 
tential 



1 



H 



5> + E ^ + E 



y ijk 



(2.1) 



i<j<k 



Since our calculation is carried out in configuration space, 
we need the nuclear interaction projected out in an opera- 
tor format allowing for calculation of two- and three-body 
matrix elements. We chose to use the operator format of 
the Argonne and Urbana family potentials. The Argonnc 
vl8 model is one of a new class of NN potentials that 
accurately fit both pp and run scattering data up to 350 
MeV with a % 2 /datum near one. The Urbana IX po- 
tential M includes a long-range two-pion exchange and a 
short-range phenomenological component. The strength 
constant of the short-range component is adjusted to re- 
produce the binding energy of the three-nucleon system. 

In configuration space, the many-body correlated 
ground state, |0), of H is written as a linear combina- 
tion of all possible n particle configurations, which are 
defined in terms of the npnh creation operators, (Oq = 1, 
0{ = a^a/^, 0\ = aj^a^a^a/jj acting on the refer- 
ence state, |0). Since the nucleon system obeys the Fermi 
statistics, the natural choice for the vacuum |0) is a n 
single-particle Slater determinant. In the coupled-cluster 
expansion, we have 

|0) = e s1 |0), 

where is the cluster correlation operator, defined in 
terms of its p/i-creation operators expansion as 
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In order to determine the amplitudes S n we solve a set 
of non-linear equations, which may be obtained using a 
variational principle 0]: 
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n > 1 



(2.2) 



The coefficients of the operators {S„} carry the physical 
significance of nuclear correlations. 

In reference (3| we have outlined our approach to solv- 
ing the coupled-cluster equations (2.2). We first intro- 
duce a mean-field Hamiltonian in terms of the one-body 
kinetic energy operator and a mean-field potential. The 
mean-field is partially constrained by the assumption 
that Si = 0, or the maximum overlap basis. The set 
of single-particle energies and wave functions which de- 
termine the configuration space is determined by diago- 
nalizing the mean-field Hamiltonian. 

The mean-field component calculation is carried out 
together with the calculation of the 2p2h correlations 
in a self-consistent manner. Even though the 3p3/i and 
4p4h correlations are not calculated explicitly, the 2p2h 
correlations are implicitly corrected for the presence of 
higher-order correlations, consistently with our trunca- 
tion scheme. In this respect we replace the traditional 
truncation in S„ Q by a truncation in l/e„, where e n is 
the npn/i excitation energy. The results presented here 
are obtained by keeping only terms up to second order 
in 1/e. 



III. CENTER-OF-MASS HAMILTONIAN 

The correlated ground-state |0) is a function of A 
sets of coordinates {fk,k — 1 . . .A}. Consequently, the 
ground state wave function is not translationally invari- 
ant. As a result, it is common practice |^] to attempt 
solving the Schrodinger equation for the ground state of 
the internal Hamiltonian 
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instead of the Hamiltonian (2.1). The internal Hamil- 
tonian is entirely written in the center-of-mass frame 
by removing the center-of-mass kinetic energy, Tqm = 
Pq M /(2iti A), with m being the nucleon mass. Both the 
two- and three-nucleon interactions are given in terms of 
the relative distances between nucleons, so in this respect 
no corrections are needed. 

However, this does not prevent the contamination of 
the correlated ground-state |0) by the center-of-mass mo- 
tion. Especially since we are interested in the calculation 
of the excited states spectrum, this contamination will re- 
sult in spurious excitations of the center-of-mass. More- 
over, the calculation of observables to be compared with 
experimental data which are taken in the center-of-mass 
frame, requires the evaluation of a many-body expansion 
as explained in Q . This in turn relies on the assumption 
that we can neglect the correlations between the inter- 
nal and center-of-mass degrees of freedom at the level of 
the correlated ground-state |0). In order to ensure such 
a separation, a supplemental center-of-mass Hamiltonian 
is added 
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The center-of-mass Hamiltonian has the role of constrain- 
ing the center-of-mass component of the ground state 
wave function Then, we have to identify a domain 
of values for the parameters 0cm and fl 2 such that the 
binding energy 
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is insensitive to the choice of these parameters. 

For the purpose of this analysis we leave out the three- 
nucleon interaction, and show results when including 
only the vl8 NN interaction. Figures |l| and || show the 
Pcm dependence of center-of-mass and binding energies, 
respectively. 



IV. THREE-NUCLEON INTERACTION 
CORRECTIONS 

We will present here the corrections necessary to take 
into account the three-nucleon interaction as part of the 
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general nuclear interaction. As discussed before we con- 
sider the three-nucleon interaction as the sum of two com- 
ponents: a long-range two-pion exchange and a short- 
range phenomenological component. The calculation of 
the three-nucleon interaction has been presented else- 
where [Hoi. 



Given the form (2.1) of the Hamiltonian, we write the 
operator Vtni m second quantization as 
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Here, the matrix elements are given as integrals involving 
the single particle states (including spins) 

Va 1 b 1 c 1 ,a2b2C2 

= ^ Ol (l)0 6l (2)0 Cl (3)|F tni |^ 2 (l)^ 2 (2)0 C2 (3)). (4.2) 

In addition to being symmetric with respect to the inter- 
change of the particles labelled 2 and 3 



' aibici ,a2&2C2 *ai ci bi ,a2C2&2 



the integrals (4.2) also satisfy the symmetry 
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The last property is a consequence of the cyclic sums 
involved in the definition of the interaction, which make 
the interaction invariant with respect to the labelling of 
the particles. 

In previous applications the approximation has been 
made that such an interaction can be represented by a 
density dependent two-body interaction. While such a 
substitution is the easiest modification, it has been stated 
that this is insufficient R|. However, the rigorous in- 
clusion of the three-nucleon interaction in configuration 
space using the coupled-cluster method is seriously ham- 
pered by our present computing capabilities and the nec- 
essary size of the configuration space. I dea lly, we would 
like to calculate all integrals of the form ( |4.2| ) without any 
artificial restrictions. In practice though, we must limit 
ourselves to calculating matrix elements of the form 
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where a and f3 couple to spin only. We are then faced 
with a compromise: Since matrix elements of the form 
(4.5) are all we need in order to calculate exactly the 
first- and second-order (S2) contributions to the mean- 
field and binding energy, the leading orders in our ex- 
pansion are treated rigorously correct. Then we make a 
reduction of the three-nucleon interaction to an effective 
two-body interaction, and use this effective interaction 
when dealing with higher-order corrections (S ra , n >3). 
This is achieved by defining the effective two-body inter- 
action as 



V 



tni.den . 



tni.den 



4 V a 1 b 1 .a 2 b 2 a ai a b 1 a b2 a a2 I l 4 - D J 

aifcia 2 fe2 



where we define the matrix element V„ K „ h to be equal 

a\b\.d2b2 ^ 



to 
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The definition (4.7) has been inspired by the form of 
second order (S2) contributions to the binding-energy 
and has the additional advantage of being fully anti- 
symmetric, so that all procedures developed when dealing 
with the nucleon-nucleon interaction [ pd| , can be natu- 
rally extended to handle the three-nucleon interaction. 
Note that the first two terms in Eq. ( [4.7[ ) are equiva- 
lent to the standard density-dependent reduction of the 
three-body force. 

We shall now detail the changes necessary to take into 
account the effects of the three-nucleon interaction in the 
calculation of the binding-energy and mean field using 
the coupled-cluster formalism. 



A. Binding Energy Corrections. 

We are only interested in the total binding energy when 
the wave function satisfies the Hartree-Fock conditions. 
Thus, it suffices to compute 

(E) = (0|T|0) + (0|V|0> + (0|S 2 V|0) + <0|S 3 V tni |0) . 

(4.8) 

The first order corrections to the binding energy are due 
to the expectation value of the three-nucleon interaction 
in the uncorrelated ground state. We have 

(0|V tm |0) = - I Vhih 2 h 3 , h!h 2 h 3 - Vh 1 h 2 h 3 , hih 3 h 2 

~t~ *h\h2h 3 . h 3 h\h2 *hi hshs ,h 3 li2 hi 

+ Vh 1 h 2 h 3 Ji2h 3 h 1 - Vh 1 h2h 3 ,h 2 hih 3 \ ■ (4-9) 



Using the symmetries ( |4.3] , 4.4), the last equation be- 
comes 
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(4.10) 
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We notice that on ly the first two terms in Eq. ( 4.10| ) 
have the form ( |4.5| ). However, in this particular case, it 
is a simple endeavor to calculate the missing third matrix 
element. We find that the magnitud e of this term is small 
compared to the sum of the terms (4.5). 

Second order contributions are calculated exactly as 

(0|S 2 V|0) = \ £ S PlhuP2h2 V^ hi . (4.11) 

Pih 1 p 2 h 2 

At the present time, the third order corrections have not 
been evaluated. Their magnitude is expected to be small, 
and their inclusion is definitely next-order in terms of the 
present expansion. 



B. Mean Field Corrections. 



of the three-body force, Eq. ( |4.7| ), for the account of the 
3p3/i correlations. Then, by making use of the full anti- 
symmetry of S2 and 'Vtni.den we can show that the re- 
quired corrections can be written as 



PlPlh\h 2 \ 
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To maintain the symmetry we have also added simi- 
lar terms to the hh- and pp-terms of the single particle 
Hamiltonian. 



In our approach to the coupled cluster formalism, the 
single particle orbits are cigenfunctions of an mean-field 
Hamiltonian, defined as the sum of a one-body kinetic 
energy term and a one-body mean-field potential. The 
later is not unique, and in the maximum overlap hypoth- 
esis Si = 0, the mean-field is defined as 

(0|U |lplfe> = (0||v+ [S a ,V] + [S 3 ,V]}|lplfc>. 

(4.12) 

Correspondingly, the contributions due to the three- 
nucleon interaction can be written as the sum of three 
terms 

(0|Vi m at a/l |0) 

+ (0|[S 2 ,V tm ]at a „|0) + (0|[S 3 ,V tni ]at a/t |0). (4.13) 

In leading order, the three-nucleon interaction correction 
of the mean-field is given as 

(0|V tni aj,a/i|0) = - ^2 \Vh 1 h 2 h-h 1 h 2 p ~ Vh 2 hh 1 ,h 2 h 1 p 
h x h 2 

~\~'hih 2 ,h' 1 h 2i p 1 hi *hi h 2 h,h± p h 2 
~t~ ^ / /i2 h h\,h\h 2 p *hih 2 h,h 2 h\ p 



Again, using the symmetries (4.3, 4.4), Eq. (4.14) be- 
comes 

(0|Vf„j ata^O) = 1 2 Vhihz h,h 1 h 2 p — Vh 1 hh 2 Ji 1 h 2 p 



~t~ vhih 2 h,h 2 p hi r> *hi h 2 h,h 2 h\ p 
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V. OBSERVABLES 

The results we report here are carried out in a 50hu> 
configuration space. In general, the expectation value of 
any arbitrary operator A in the correlated ground-state 
|0) is obtained as 

a= (0|e s Ae- s S+ | 0) , 

where similarly to the cluster correlation operator , the 
new operator £>t is expanded out in terms of p/i-creation 
operators 

n 

The amplitudes S n are calculated iteratively in terms of 
the S n amplitudes ||. 

The resulting binding energy and charge radii for the 
Argonne ul8 with/without the Urbana IX potential, are 
shown in Table |. Note that the charge radii correspond 
to the "theoretical" charge density as explained below. 
The calculation is shown to be stable with respect to the 
two cut-off parameters, l ma x and n maxi which control the 
size of the configuration space as N max = 2n max + l max . 
Figure || shows the binding energy for 16 O calculated 
using the vl8 and Urbana IX interactions, as a function 



(4.14) of the (Z 



max j IL max 



) cut-off pair, and show the calculation 



The second- and third-order contributions in Eq. ( 4.13| ) 
look very similar when one uses the proposed reduction 



to be reasonably converged for l max — 11, 12 and n Ti 
22, 25, same as the wl8 based calculation of jj). 



A. Charge Form Factor 

The charge form factor is calculated as the sum of 
a one- and two-body operator. The one-body compo- 
nent is evaluated in the nonrelativistic one-body Born- 
approximation picture, where the charge form factor is 
given as 
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Here fk(q 2 ) is the nucleon form factor, which takes into 
account the finite size of the nucleon k. In our calcula- 
tion we use the Iachello-Jackson-Landc nucleon-form 
factors. Note however that at the rather small values of q 
involved in the present discussion, q < 4 fin , the dif- 
ferences between the various models of the nucleon form 
factor are expected to be small |u|. 

As seen from previous calculations Q, the main con- 
tributions from the two-body charge density not related 
to center-of-mass corrections is expected to come from 
the the tt- and p-exchange "seagull" diagrams. These are 
taken into account as referred in the model-independent 
part of the "Helsinki meson-exchange model" jl3| . In this 
model the pion and p-meson propagators are replaced by 
the Fourier transforms of the isospin dependent spin-spin 
and tensor components of the vl8 NN interaction, in or- 
der to ensure that the exchange current operator does 
satisfy the continuity equation together with the interac- 
tion model. As seen from Fig. |], the contributions of tt- 
and /3-exchange give a measurable correction for q > 2 
fm _1 . The order of magnitude of the n- and p-exchange 
contributions compares well with the calculation of ||] 
for 16 using the Argonne wl4 potential and Urbana VII 
three-nucleon interaction. 

Since the correlated ground state wave function |0) is 
not translationally invariant, extra care needs to be taken 
in order to account for the effect of the center-of-mass 
motion on the expectation value of the operator in the 
ground state. Center-of-mass corrections have been dis- 
cussed in H , and rely on the assumption that the motion 
of the intrinsic coordinates f" k — — R cm , k = 1 ... A, 
and the center-of-mass are not correlated, and that our 
correlated ground-state |0) provides indeed a good de- 
scription of the internal structure of the nucleus. 

Figure [| depicts the center-of-mass corrected charge 
form factor both in the impulse approximation and in- 
cluding the meson-exchange two-body component, to- 
gether with the experimental results for q < 4 fm -1 [ fli) . 
These various approximations of the theoretical charge 
form factor correspond to the v18+XJIX.+Hcm interac- 
tion in Table Q. 

Figure [5] shows the corresponding charge density to- 
gether with the "experimental" charge density as derived 
in p4] |. This calculation simply represents the result of 
Fourier transforming the charge form-factor presented in 
Fig. ^. However, in doing so we use predictions for a mo- 
mentum transfer greater than 4 fm _1 , where the present 
calculation is not expected to offer reliable predictions. 
This results in uncertainties of the charge density distri- 
bution at short distances r. 

Also, the distortion effect related to the interaction of 
the electron probe with the nuclear Coulomb field has 
not yet been accounted for. Thus, we use the Distorted 
Wave Born Approximation (DWBA) to calculated a new 
charge form factor Q , which is depicted in Fig. ||. Note 



this calculation is now sensitive to the energy of the inci- 
dent electron as shown in Fig. ^. This time the agreement 
is quite impressive, and the calculation of the correspond- 
ing charge density will give a result considerably closer to 
the "experimental" one. As a result of this exercise, we 
conclude that it is better to compare the results of the 
distorted charge form factor calculation with the experi- 
mental charge form factor over the domain of momentum 
transfer covered by the experiment. In any event the 
charge form factor represents the primary result of the 
experiment, and comparing form factors is thus more re- 
liable than comparing model- dependent charge densities. 

Note that the persistent discrepancy at large q may 
be a signature of the break down of the nonrelativistic 
approach to calculating the charge form factor. However, 
due to the various other approximations made so far, a 
calculation based on the relativistic description of the 
one-body charge form factor E^ j may be too expensive 
at this time. 



B. Coulomb Sum Rule 

The Coulomb sum rule, Sl(q), represents the total in- 
tegrated strength of the longitudinal response function 
measured in inclusive electron scattering. The Coulomb 
sum rule is related to the Fourier transform of the 
proton-proton distribution function in the nuclear ground 
state [jl7| . As such, this quantity is sensitive to the short- 
range correlations induced by the repulsive core of the 
NN interaction. In the nonrelativistic limit we have E^] 

SL(q) = ^(6\p i p(q)\0)-^\(0\p(q)\d)\ 2 



= l + p LL (q)-^\(0\p(q)\0)\ 2 , 
where p(q) is the nuclear charge operator 
1 A 

p(q) = 2 E e ^ (1 + ^.0. 



(5.1) 



(5.2) 



In Eq. (5.1), the longitudinal-longitudinal distribution 
function Pll(<z) is given in terms of the proton-proton 
two-body density as 

PLL(q) = dn dr 2 j (q\ri - r 2 \) p {p ' p) (fi, r 2 ) ■ 



(5.3) 

We have outlined recently ||] the calculation of the 
proton-proton two-body density in the coupled-cluster 
framework, as the expectation value 

P {p ' p) (h,r 2 ) (5.4) 

= \ E( 8 l^l - Wfa - + T z ,i)(l + T z ,j)\0) , 
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with the normalization 



dn I dr 2 />' p )(fi,f 2 ) 



(5.5) 



Figure ^ depicts the results of the present calculation 
for 16 0. Agai, the theoretical calculation corresponds to 
the v18+\JIX+Hcm interaction in Table |. Since no ex- 
perimental data are available for 16 O we compare these 
results with the 12 C experimental data of |jl9), where 
an estimate j20| for contributions from large u> has been 
added. Preliminary theoretical results for 12 C obtained 
using the coupled-cluster method have been recently re- 
ported in [^l| and shown to be close to the 16 O result. 
The large error bars on the experimental data are largely 
due to systematic uncertainties associated with tail con- 
tribution |E2fl. 



VI. CONCLUSIONS AND OUTLOOK 

The goal of our effort is to make a contribution to 
the ongoing effort of building realistic models of nu- 
clear structure that explicitly account for realistic cor- 
relations. In a first stage we focus our attention on ob- 
taining a realistic description for the ground state of a 
doubly-magic nucleus using the coupled-cluster method. 
In this sense, the present calculation represents the most 
detailed calculation available today, using the coupled- 
cluster method for a nuclear system with A > 8. We 
base our calculation on what is considered the best real- 
istic description of the two- and three-body interactions 
available today ^ . 

At this stage we have proved conclusively that it is 
possible to choose a large enough configuration space to 
handle the relatively hard-core of the nucleon-nucleon in- 
teraction. We believe that we have also identified the 
dominant contributions of the three body interaction, 
and we have introduced an effective two-body interac- 
tion, accordingly. 

An interesting conclusion from the charge form factor 
calculation has to do with the importance of the distor- 
tion effect due to the interaction between the electron 
and the Coulomb field of the nucleus. In this context we 
have shown that it is desirable to compare the distorted 
charge form factor with the experimental one, and to de- 
emphasize the comparisons regarding the nuclear charge 
density, which is not very surprising as the charge form 
factor represents the primary result of the experiment. 

In this paper we have confined ourselves to presenting a 
purely nonrelativistic description of the charge form fac- 
tor in 16 O. In the future it will be interesting to get a bet- 
ter understanding of the importance of relativistic effects 
on the charge form factor calculation. In this sense we 
plan to redo this calculation using the relativistic descrip- 
tion of the one-body charge form factor of , together 
with a consistent relativistic derivation of the two-body 
meson-exchange density [23]. 



With the calculation of the 16 O ground state completed 
we intend to extend our formulation to address the cal- 
culation of discrete excited states as well as neighboring 
odd-even nuclei. This is a necessary step in the quest of 
modeling the (e,e'N) reaction, where the final state has 
the asymptotic form of a distorted wave times a discrete 
state of the (A-l) nucleus. However, this is not a solution 
to the Hamiltonian close to the origin, and thus the wave 
function needs to be modified in the region of the origin. 
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TABLE I. Energy expectation values and theoretical 
charge radii for 16 O. 



Interaction 


B.E. 


r.m.s 




[MeV/nucleon] 


[fin" 1 ] 


vl8 


-5.9 


2.85 


vlS + UIX 


-7.7 


2.74 


vl8 + UIX + Hcm 


-8.0 


2.62 


expt. 


8.00 


2.73 ± 0.025 
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FIG. 4. Charge form factors for le O, using the Argonne 
vl8 and Urbana IX potentials. 
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FIG. 2. Pcm dependence of binding energy (using vl8). 
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FIG. 5. Charge density for 16 0, including meson exchange 
contributions. 
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FIG. 3. (Imax ,rimax) cut-off dependence of binding energy 
(using vl8 and UIX). 
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FIG. 6. Charge form factor for ie O, obtained in the DWBA 
picture («18 + UIX + ME). 
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FIG. 7. Coulomb sum for 1(i O, compared with "exper- 
imental" 12 C data which include theoretically determined 
high-energy tail corrections. 
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